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Abstract 

In this paper, we study the vanishing viscosity hmit for a coupled Navier-Stokes/Allen-Cahn 
system in a bounded domain. We first show the local existence of smooth solutions of the 
Euler/ Allen- Cahn equations by modified Galerkin method. Then using the boundary layer 
function to deal with the mismatch of the boundary conditions between Navier-Stokes and 
Euler equations, and assuming that the energy dissipation for Navier-Stokes equation in the 
boundary layer goes to zero as the viscosity tends to zero, we prove that the solutions of the 
Navier-Stokes/Allen-Cahn system converge to that of the Euler/Allen-Cahn system in a proper 
small time interval. In addition, for strong solutions of the Navier-Stokes/Allen-Cahn system in 
2D, the convergence rate is cv^^"^. 
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1. Introduction 

In this paper, we are concerned with the vanishing viscosity limit for the following Navier- 
Stokes/Allen-Cahn system in x (0, -|-oo): 

Ui (u ■ V)u + VPi = i/Au - AV • (V0 (g) V0), (1.1) 

V-u = 0, (1.2) 

0i + (u.V)</. = 7(A0 -/(,/.)), (1.3) 

with initial data 

u(x, 0) = uo(x), (/)(x, 0) = (/)o(a:), x G 0, (1.4) 

and boundary conditions: 

u(x,t) = 0, dn<p{x,t)=Q, (x,t) G X (0,+oo). (1.5) 
Here, C in = 2, 3) is a bounded domain with smooth boundary 50. 
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System ()l.ip - (|1.5p can be viewed as a phase field model, which describes the motion of a 
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mixture of two incompressible viscous fluids with the same density and viscosity (see 
The two fluids are macroscopically immiscible and separated by a thin interface, (jl.ip is the 
linear momentum equation, where u is the velocity field of the mixture, (j) ^-nd Pi denote the 
phase function and pressure, respectively. V0(8> V</> denotes the induced elastic stress, which is a 
71 X n matrix whose (i, j)-th entry is §^-§^ for 1 < i,j < n. ()1.2p implies the incompressibility of 
both fluids in the mixture and ()1.3p is the phase equation of Allen-Cahn type. f{4>) = ^■(0'^ — </>). 
z^, A, 7, e are positive constants, representing the kinematic viscosity, the surface tension, the 
mobility and the width of the interface, respectively. u\qq = is non-slip boundary condition. 
And dn(l)\dn = means that the diffused interface perpendicularly contacts the boundary of the 
domain, where n is the outward unit normal to the boundary dQ. 

From another point of view, system ()l.ip - (jl.3p is closely related to liquid crystal model, 
Magnetohydrodynamics (MHD) equations, and viscoelastic system with infinite Weissenberg 
number, see 

System (|l.ip - (|1.5p has a basic energy law 

~ (\u\' + + _ i)[2\ ^ _(^[vu|2 + A7|A0 - (1.6) 
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where | • | denotes the L^.^orm. (|1.6p can be derived by multiplying (jl.ip by u, ()1.3p by 
A(—A</) + /((/))), then adding them up and integrating over Cl, also using ()1.2p and (|1.5p . Thanks 
to this important energy law (|1.6p . it can be proved that global weak solutions of system (jl.ip - 
(jl.Sp exist in both 2D and 3D case. With regular initial data, strong solutions exist globally in 
2D, but locally in 3D in a short time interval [0,T^), where T^, depends on v. Hence we just 
consider weak solutions of ()l.ip - (|1.5p in the 3D case in the proof of vanishing viscosity limit. 

At the limit case, namely = 0, Navier- Stokes/ Allen-Cahn system formally becomes the 
following Euler/ Allen-Cahn system: 

vj + (vV)v + VP2 = -AV-(VV'«)VV'), (1.7) 
V-v = 0, (1.8) 
^^ + (v.V)V = 7(A^-/(V)), (1.9) 

in X (0, +oo) with initial data 

v(x,0) = vo(x), '4>{x,0) = ijjQ^x), X G r2, (1.10) 

and boundary conditions: 

v-n(x,t) = 0, dn'4^{x,t) =0, {x,t) e dn X {0,+oo). (1.11) 

For simplicity, we take vo(x) = uo(x) and ipo{x) = <l)o{x). By taking the limit of weak formula- 
tion of the problem (1.1)-(1.5) as z^ — )• 0, we obtain the weak formulation of (1.7)-(1.11), from 
which we can derive the system (1.7)-(1.11). Thus the boundary condition of v is v • n\g^ = 0, 
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which describes that the boundary is impermeable. The basic energy law of system ()1.7p - (|l.lip 
is 

The main purpose of this article is to show that the solution of the viscid system ()l.ip - ()1.5p 
converges to that of the inviscid system (|1.7p - p.lip in a short time period as the viscosity goes 
to zero. The inviscid limit helps us to understand turbulent phenomena governed by the viscous 
equations, see Navier-Stokes equations for an example (see e.g. (ll. ItI. IsI. IiqI. S. [2^). 

Let us first recall the classical issue of vanishing viscosity limit of Navier-Stokes equations 
in a bounded domain. In order to deal with the mismatch of the boundary conditions between 
Navier-Stokes system and Euler system, Kato in introduced a boundary layer function defined 
on a boundary strip Tcu with width cv, and proved that the solution u'^ of Navier-Stokes equation 
converges to that of Euler equation in L^-norm, provided 

vj |Vu^(t)||2(r^^)dt ^ 0, as 1/ ^ 0, (1.13) 

where [0, T] is the interval on which the smooth solution of Euler equation exists. Later, Temam 
and Wang in [23] and 25[ improved Kato's condition by replacing the total gradient in ()1.13p 
with tangential derivatives only, but required a slightly thicker boundary layer. 

With the same boundary layer function inl?^, and under the same assumption ()1.13p . we 
can also prove that weak solution of ()l.ip - ()1.5p converges to smooth solution of (|1.7p - ()l.lip in 
the basic energy space in [0, T] as the viscosity tends to zero. In particular, for 2D case, strong 
solution of ()l.ip - ()1.5p converges to smooth solution of ()1.7p - ()l.lip at the rate cz/^/^. The main 
difficulty in the proof is to deal with coupled terms of u, v, and ip. A key point is that the 
highest order coupled terms can be canceled, similar to the derivation of the energy law (|1.6p . To 
estimate other coupled terms, high regularity of solutions to ()1.7p - ()l.lip is required. Therefore 
we first establish the local existence and uniqueness of smooth solutions to Euler / Allen-Cahn 
system ()1.7p - ()l.lip . This proof requires the energy law ()1.12p . so we employ the modified 
Galerkin method introduced in 11[, namely, we only project (jl.7p into finite dimensional space. 



but solve p.yp for any given v by the fixed point theorem, such that the approximate solutions 
still satisfies the energy law ()1.12p . 

This article is organized as follows. In section 2, we introduce mathematical preliminaries 
and state our main results. In section 3, we establish the local existence and uniqueness of 
smooth solutions to the Euler/Allen-Cahn system. In section 4, we prove that the solution of 
the Navier-Stokes/ Allen-Cahn equations converges to that of the Euler/Allen-Cahn equations 
under the condition ()1.13p . In the last section, we show some results for related systems. 



3 



2. Preliminaries and Main Results 



Let C^j^(f]) be the space of all divergence free vectors in (C^(17))" (n = 2,3). We denote 
by H the closure of C'^^^(O) in (L^(r2))". Moreover, we set 

V = {H^{n)r riH, V2 = iH\n)r n v, 

W = {(pe C°°(J1), dn(l)\dn = 0} , ^s = the closure ommH'{n), s G N+, 

O = {v € {C^{n))'', V • V = in V n = on af)} , 

Xo = the closure of O in (l2(J7))", = n (//"(l)))", s e N+. 

Obviously, C Xq (continuous embedding), Xg CC Xq (compact embedding) for s > 1. 

Let be the boundary strip of width 5. We shall take 6 = cu with c > being a small 
constant. We denote | • | and (•,•) as the norm and scalar product in L^{Q) or (L^(0))". For 
any positive integer s, we take | • \s and ((•, ■))s as the norm and scalar product in H^{0,) or 
(i7*(r2))", {{f,g))s = Yl\a\<s {^"'fiD'^d) where D° is a multi-index derivation. We also denote 
I • |x as the norm in other Banach spaces X. 

Following almost the same arguments as in 11[, for any fixed > 0, using the important 
energy law (|1.6|) . we have the following well-posedness result for the Navier-Stokes/Allen-Cahn 
system: 

Theorem 2.1. Assume that n = 2,3. uq G H and (/)o € $1- Then the system (|l.ip - ()1.5p has a 
global weak solution (u, (j)) such that for all T £ (0, +00) 

u G L^{0,T;H)nL'^{0,T;V), 
(/)EL~(0,r;$i)nL2(o,r;$2). 

Moreover, if uq £ V and (po G ^2, then in the 2D case, system ()l.ip - ()1.5p admits a unique global 
strong solution such that for all T G (0, +00) 

ue L'^{0,T;V)nL^{0,T;V2), (2.14) 

<^GL~(o,r;$2)nL2(o,r;$3); (2.15) 

while in the 3D case, there exists aT^ > depending on |uo|i, |</'o|2 CLiT'd v such that (|l.ip - (|1.5p 
has a unique strong solution in [0, Tjy) and ()2.14p - (|2.15p hold for T < T^. 

For weak solution (u, 0) to system (|l.ip - (|1.5p . we easily deduce from energy law (|1.6p that 
for any T > 0, 

V [ |Vu[2(t)dt <C, [ \(l)\lit)dt < C, (2.16) 
Jo Jo 

where C depends on uo,(/)o and 

We also prove the local existence of smooth solutions to Euler/Allen-Cahn system. 
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Theorem 2.2. Assume that n = 2,3. vq € Xg, ipQ € ^s+i, where s > 3 is an integer. Then 
there exists a small constant > depending on |vo|s and |^/)o|s+i; such that the system (jl.7p - 
(jl.lip admits a smooth solution (v, P2)'0) [O,^!,,] satisfying 

veL°°(0,r,;X,), P2eL'^{0,n;H'+\n)), V e i~(0, T,; (2.17) 

Moreover, (v^Tp) is unique in the sense o/ (L°°(0, T^,; L^(il)), L°°(0, T^,; //^(il))) -norm. 

Finally, we show the convergence of solutions of Navier-Stokes/Allen-Cahn system to that 
of Euler/Allen-Cahn system. 

Theorem 2.3. Assume that n = 2,3. Let {u,(j)) be a global weak solution to ()l.ip - ()1.5p . and 
{v,ip) a local smooth solution to (jl.7p - (jl.lip on 0, x [0, T*]. Assume that 

vj |Vu|^2(r^_^)(t)di ^ as ^ 0. (2.18) 

Then 

|u(i)-v(t)p + A|(/)(t)-?/;(t)|?^0 /or a.e. t G [0, T*]. (2.19) 

Moreover, if n = 2 and (u, (/>) is a global strong solution to ()l.ip - ()1.5p . then the assumption 
(j2.18p is automatically satisfied, and the convergence rate in (j2.19p is cu^^'^ . 

Remark 2.1. /n general, we do not know whether the weak solution {u,(f)) satisfies the assump- 
tion ()2.18p . We only know Jq* |Vu|^2(p ^{t)dt is bounded for weak solutions (see ()2.16p ). But 
in 2D, the global strong solutions exist. If{u, (f)) is a global strong solutions, i.e. || Vu||2,cx)(o,T;L2(q)) 
is bounded and its upper bound is independent of v, then u satisfies (|2.18p obviously. 

Remark 2.2. The above theorems remain true if we replace the Neumann boundary condition 
for (p and ip with homogeneous Dirichlet boundary condition. 



3. Local smooth solutions to Euler/Allen-Cahn system 

In this section, we give the proof of Theorem 12.21 in the 3D case. The 2D case can be proved 
similarly. We shall apply a modified Galerkin method, followin g th e spirit of 1]|. We choose 
an orthonormal complete basis {^k}'k'=i C (s > 3) used in [20(] for Euler equation, which 
satisfies 

f ((wfc, v)), = Afe(wfc, v), Vv G Xs. 
\ {■Wk,MVj) = 6kj,i.e.,{{-Wk,Wj))s = XkSkj- 

Here, < Ai < A2 < • • • < A„ < • • • and A„ — 00 as n — )• 00. 

For fixed integer m > 0, denote Xsm as the space spanned by {wi,W2,--- ,Wm}- Let 
Pm ■ Xq — 7> Xsm be the orthonormal projection. We consider the approximate problem: 

= -Pm (-(Vm • V)Vm - AV • (V^/^m (g) VV'm)) , G ^sm, (3.21) 
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^ + (v^ • V)^^ = 7 (AV'm - fi^Pm)) , (3.22) 
Vm(0,x) = PmVo(x), Vm(0,x) = ^0(2;), (3.23) 
Vrn ■ nlsQ = 0, dn'ipm\an = 0. (3.24) 

Here, we only project ()1.7p into finite dimensional space Xsm, but do not project ()1.9p . We 
shall see later in ()3.28p that, the approximate solutions {vnni^m) still satisfies energy law. This 
guarantees the global existence of {vm,ipm)- 

We apply Schauder fixed point theorem tp_deduce the local existence of solutions to ()3.2ip - 
()3.24p . The arguments are similar to that in (ill . Therem 2.1]. For readers' convenience, we give 
the outline of the proof here. Consider a closed, convex subset in (C[0,To])™: 



1 -,2 ... mN 



^ m \ 1/2 

l5m(^)P < M, < t < To; gj„(0) = (vq, w,), i = 1, 2, • • • , m. 



where M and Tq are positive constants to be determined later. Given a {gj^,g'^, ■ ■ ■ ,5™) G D, 
we get a Vm(x,t) = EI^i 5'm(*)wi(x) satisfying |vm,|ioo(f^) < c^M for t G [0,ro], where Cm 
is a constant depending on m. With this in ()3.22p . there exists a regular solution iprn on 
n X [0,To], and 

|VV™P(t) <e'=^^'*(^|VV^o|' + 2^F(V'o)d^) =:e^^*^'*C7(lV^o|i,|J^|), [CTq]. 

Then we substitute this into ()3.2ip and look for a solution = X^ilLi 5m(*)'^i(^)- I* is 
standard that (|3.2ip is equivalent to an ODE system of glnit) (* = 1) 2, • • • , m), with initial data 
5rn(0) = (vo,Wj). There exists a local solution {gl^,g^,--- ,5^) ^ (C'^IO, To])*". Moreover, it 
turns out that 

= (EI(vo,w.)p)''% ^^^"^;J/;'^'"^ (e---^^^-- - 1) (3.25) 

for t € [0, To]. Let M = 2 i^Jl^i |(vo, Wj)p)^''^ + 2. Then there exits a small Tq such that the 
right hand side of (13:251) is not more than M. Therefore, {Ya=i \9mif')?f''^ < M for t G [0,To]. 
Thus {gl^,g^,--- ,g™) G D. And the mapping £ : (c/^,c/^,-- - ,5™) h^. {g}^,gl^,--- ,gj^) is a 
compact operator in (C[0,Tq])™'. Applying Schauder fixed point theorem, C has a fixed point. 
Thus we obtain the following lemma: 

Lemma 3.1. There exists qTq > 0, depending on vojipQ^m and 0,, such that (|3.2ip - (|3.24p has 

a weak solution {vm,4'm) in 0, x [0,Tq]. 

In order to show the global existence of weak solution to ()3.2ip - ()3.24p . we establish a priori 
estimates next. Assume (v^, ipm) is a weak solution to (j3.2ip - (j3.24p in 17 x [0, T] for certain T > 
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0. Taking the inner product of ()3.2ip with in L^(i7), using the fact that ((v-m • V)vm, v.^) = 
and V • (VV'm Vipm) = V (^\^^^ + A^'^VV'm, we get 

^^|VmP = -A (AV'mV^/^m, Vm) . (3.26) 

Taking the inner product of p.22p with 7 {Aipm — /(V'm)) ™ -^^(^)) we have 

4(^|VVm!^ + A/ +A7|A^„-/(V„^)P = A((v^•V)V'™,AV'™), (3.27) 



where we used {{vm ■ y)'ipm, f{ipm)) = fn^rn ■ VF(V'm) = where F{ipm) = -^{tpm - 1)^- 
Adding the above two resultant, we obtain 

4 f ilvmP + ^IVV-^P + A / F{i;m)dx] + Xj\Ai;m - /(V-m)!' = 0. (3.28) 



dt \2 

Hence, 



sup (|v^|2(t) + A|VV'mP(t) + A / F{^Prr,)dx] < IvoP + AlVV'oP + A / F{iPo)dx. (3.29) 
o<t<T \ Jn / Jn 

By Lemma l3.ll and the above inequahty, we have 

Theorem 3.1. For any integer m > 0, vq € Xq, ipQ G $1, system (|3.21|) - (|3.24|) has a weak 
solution {vm,ipm) inVtx [0, 00). 

In what follows, we do higher order energy estimates, so that we can pass to limits as m — > 00. 
Firstly, sincG G Xgyyi^ we can rewrite fj3.2ip cis 

dw 

= -(Vm ■ V)Vm - \A'4)m^llJm " Vg^, (3.30) 

Ot 

for a certain qm- By taking divergence operator on both sides of (j3.3U|) on 0, and taking scalar 
product of both sides of (|3.30p with n on dVt, we know that qm satisfies 

Aqm = -AV • {A^pm'^tprn) " ^ VjVmjVjVmj, in Q 

— — = -XAipm^ \- } ^rmi^^rnj^ij, on dn 

on on ^-^ 

where the function Tij depends only on dQ (see {20! . Lemma 1.1]). By elliptic estimates (see also 
20, Lemma 1.2]), for s > |, 



lVg,„|, < C (lAV^^UWml. + Wml'i) ■ (3.31) 

Taking inner product of p.30p with XkOkm^k in L'^i^)^ SLnd adding in A;, A; = 1, 2 • • • , m, we 
obtain 

2^|Vmls = ((-(Vm • V)Vm - AA^mVV'm - V^m, Vm))s, (3.32) 
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For s > 3, we estimate the righthand side of (|3.32p to get 



1 d 
2di 



Vmis < |((Vm • VVm, Vm))s| + IA((AV'mVV'm, Vm))s| + |((VQ'm, Vm))s| 

< C\Vmfs + CXlAlpmlsl'^lpmlsWmls + {"^qmlsl^mls 

< C\V^\1 + ^\AlPm\'i + ClVlPmlt + C\Vm\t + C ( | AV'rn U |, + |v„|2) |v^ 



< ^\Ai;m\l + C\V^l^m\t + C\^^m\t + C, 

where we used ((v^ • V)V*Vm, V**Vm) = and (|3.3ip . 



(3.33) 



Remark 3.1. In the 2D case, H^+%n) C L°°{n), Ve > 0. |((v„ • Vvm,Vm))s\ in ffOH]) can 

be controlled by |vmj" optimally, where 2 < a < 3. This term makes it difficult to obtain global 
existence of smooth solutions. 

Applying (i = 1, 2, • • • , s; s > 3) to ([3:22]) (here, V* = V(-A)(*-y/2 if i is odd, V* = 
(— A)*/^ if i is even), then taking the inner product of the resulting equation with W^Ail^m, 
and using the boundary condition ^^^\dn = and higher order natural boundary conditions 
such as ^ gn'" ~ ^' obtain 



Xd_ 
2di 



|V'+VmP + X7\V'A^m\^ = A (V^(v^ • V)i^m, V^A^^) + A7 (V7(V'm), V*AV„^) . 



Adding the above resultant in i (i = 1, 2, • • • , s; s > 3), we get 
A d 

\Vipm\2+Xj\Aip„i\; = A 

1=1 i=l 

Next we estimate the righthand side terms. 



|V^„|^ + A7|AV',„|^ = A^ (V'(v„ • V)V'm, V^AV'™)+A7^ (V7(V'm), V'AV'rn) . (3.34) 



A 



^(VXv™•V)V„^,V*AV„^) 
i=l 

s 

< ^A(|VV^||VVm|L- + |Vr„U-|V^+Vm!) [V^A^^^ 
1=1 

s 

< J]A(|VV^||VV'™|2 + |v™|2|V'+Vm|) IV'A^^I 



< 



1=1 



E 

i=l 

At 



^7, 



m 



< -^lAtPmli + cMt + civ^iJmfs, 



(3.35) 



At (v7(V'm),v*A^„) 

^(V^(V'^),V^AV'„) 
12At 



^(vVm,V^A^^) 



-^{^'i^m,V'Aij^) 
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= : /1 + /2+/3, 



(3.36) 



where S represents the sum with respect to integers ii, 12,^3. satisfying < ii < ^2 ^ ^3 < * — 1 
and ii + ^2 + ^3 = i- 

Using the fact that V'm is uniformly bounded in L°°{{),T\H^{^)) and C L^, we get 

< C|Vm|i6|VVm|L6|V^A^^| 

< ^IV^AV^^P + CIVV-™!?, (3.37) 



3 

< Cj][V*AV™in (C'|V^™^-"''|V*AV„^r'=+C|V^™|) 



k=l 



< C^IV^AV'™! (|V^AV„^|^ +1) , 



< ^iV'AxPml^ + C, (3.38) 
o 



where we used the Gagliardo-Nirenberg interpolation inequalities (cf. [30|): 



|V*'=Vm|L6 <C|VV™|'-'^'=|V'AV„^r'= + C|VV„^|, «fc = T^, A: = 1,2, 3. 



and ai + 02 + as = ^ < 1. 

I/3I < C|VVm||V'AV'r«| < ^IV^AV^^P + C|VVm|'. (3.39) 



Using the above estimates, we conclude from ()3.34p that 



^^IV^^I^ + ^\A^P^\l < C\V^P,n\t + C\v^\t + C. (3.40) 
Adding IK33\i and I^JOh . we conclude 

~ {Wml^s + A|VV^„|2) + ^\A^m\'i < CWmlt + C|VV^^|f + C. (3.41) 

Denote ym(i) = Wmll + AjV^'mls- (|3.41|) can be rewritten as 

jYm{t)<C,Y^{t) + C2, 



Y^O) = Iv„|2(0) + A|VV'™|^(0) < |vo|2 + A|VVo 



Hence there exists a T^, depending only on Ci,C2, and |vo|s, iV'ols+ii such that Ym(t) < N on 
[0,T^,], where > is a constant independent of m. Therefore, as m — > +00, 

Vm remains bounded in L°°(0, T*; Xg), (3.42) 
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1pm remains bounded in L°°(0, T^,; <I>s+i)). (3.43) 



By the equation ()3.2ip and ()3.22p . it is easy to know that 



dt 



remains bounded in L°°(0, T*; Xs_i), (3.44) 
remains bounded in L°°{0,T^;^s-i). (3.45) 



dt 

Using (j3.42p - (j3.45p and the standard compact embedding theorem, the passage to the hmit is 
standard. We obtain the existence of v G L°°{0,T^; Xs), ip € L°°(0, T,,,; <I>s_|_i) such that 

= (-(v-V)v-AV-(VV'®V^),^), V^eXo, a.e. t € [0,r,], 

(^,^) + ((v,n-V)V,95) = -7(VV,V(^)-(/(V'),v^), y^eH\n), a.e. t G [o,r,]. 

The proof of global existence of smooth solutions to (jl.7p - (|l.lip is finished. 

It is standard to prove the uniqueness of smooth solutions to (|1.7p - ()l.lip . We omit it here. 
The proof of Theorem 12.21 is complete. 

4. Vanishing viscosity limit 

In this section, we prove Theorem 12.31 In the proof, C is a constant independent of v. 

To deal with the mismatch between the boundary condition of u and that of v, we use 
the boundary layer function constructed in p]: 6 = div (x (f) A) where p = dist{x,dO,) is a 
distance function, 5 is a small positive constant, and x is a smooth cut-off function x ■ 
such that 

x(0) = 1, x(^) = for r > 1. 
A is a skew-symmetric matrix defined on x [0, T], such that 

div A = v on (90, A = on 50. (4.46) 

Obviously, is supported on Tg. It was proved in [2] that satisfies 

divG = in 0, = v on dQ, 

and has the following estimates: 

\e\L^<c, \e\<c6^, \e\Li<C6^, \et\<C6^, (4.47) 

|V0|lco < C(5"\ \ve\ < CS-^, \pVe\Loo < C, (4.48) 
|p^V6>|ioo < C6, \pVe\ < C6^, (4.49) 

provided that v has the regularity: 

I W[loo < C, [VvtlLoo < C, \vt\L^ < C. (4.50) 
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Suppose (u, </)) is a global weak solution to (|l.ip - ()1.5p . and {v,ip) is a local smooth solution 
to ()1.7p - ()l.lip on X [0, T*]. Denote w = u — v + = — ■;/', then w, ( satisfy the following 
equations: 

wt + (u • V)w + Vq = i^Au - (w • V)v + {6 ■ V)v + (u • V)0 + Ot 

-XACV(p - AAV-VC, (4.51) 
V-w = 0, (4.52) 
Ci + (v • V)C + (w • V)</) = jAC + {e-V)(l)-^a^^ + ^^ + ^^) + ^C, (4.53) 
w|t=o = 6'(x,0), C|i=o = 0, xen, (4.54) 
wlon(x,t) = 0, anC!9Q(2;,i) = 0, t>0, (4.55) 

where q = Pi - P2 + ||V0|2 - ||V?/^|2. 

Our goal is to estimate |wj(t) and |CI_f/i(^)) we hope they go to zero, equivalently, |u — v|(t) 
and \(f) — ip\Hi{t) converge to zero as ly tends to zero. 

Multiplying (|4.5ip by w and integrating over Q, we obtain 

|w|2(t) + i.| Vu|2 = u (Vu, V(v - 6*)) - ((w • V)v, w) + ((6> • V)v, w) + ((u • V)6', u) 

- ((u • V)e, v) + (0j, w) - A (ACV(/), w) - A (A^VC, w) . (4.56) 

Multiplying ()4.53p by AC and integrating over ^l, we get 

^||CP + AtIvcI^ + ^ X i^" + + 

= ^|Cp-A(w.V,^,C) + A(0-V0,C). (4.57) 

Notice that cp'^ + cp^ + ^p'^ = i ((/>2 + ^2 _^ (0 + ^)2) > 0, thus C {<f + H + V'^) dx > 0. 
Multiplying ()4.53p by AA^ and integrating over Q, we obtain 

^^|VCP + A7|ACP = A((v.V)C,AC) + A((w.V),^,AC)-A((0-V)(/>,AC) 

+^ (C('A^ + + AC) + ^1 VCP. (4.58) 

We observe that the term —A (ACV</>, w) in (|4.56p cancels out A ((w • V)</>, AC) in ()4.58p . which 
is originated from the energy law ()1.6p . This is very important, otherwise we cannot control the 
two terms. 

Next we estimate the other terms on the right-hand side of (|4.56p - (|4.58p . using (|4.47p - (|4.49p . 

|i/(Vu,V(v-6>))| < z.lVu|lVv| + z.|VuIi2(r,)|V6>ji2(r,) 

< ^iVup + z^l Vvp + Cu5-^ |Vu|i2(r,), (4.59) 
[((w-V)v,w)| < |Vv|loo|wP, (4.60) 
|((0-V)v,w)| < |0[|Vv[loo|w| < C5^|Vv|Loo[vi^[ 

< C|wp + 5|Vv|ioo, (4.61) 
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I((u-V)0,u 
|((u.V)0,v) 



1 







1 


0—udx 


< 


— u 


P 




P 



\pVe\ !v|ioo(r^) < C5^|Vu|i2(r^)|v|ioc(r^), 
where we used ()4.49p and the weU-known inequahty of Hardy-Littlewood (since u G Hq{Q)). 

K^t, w)| < \et\\vj\ < C6^\w\ < C|w|2 + C75, (4.62) 







1 


j (^-u ■ pV^ Ovdx 


< 


— u 






P 



|A(AV^VC,w)| < A|AV^|i4|VC|L4lw| 

^ J C7A|AV'|l4|VC|5|AC[^|w| (ifn = 2) 
\ CA|AV'|L4|VC|3|AC|t|w| (ifn = 3) 

< ^|ACp + C|VCp + C|A^|i4|w|2, 



(4.63) 



|A(wV0,C)| < A|w||V</.||C|l- 

' C|w||V(^||C|5|AC|^ (ifn = 2) 
C7|w[|V(^I|C|3|AC|^ (ifn = 3) 

^|ACP + C|w|2 + ClV0|4|CP (ifn = 2) 
^|ACP + C|w|2 + C|V(/)nC|2 (ifn = 3) 



< 



< 



(4.64) 



\xie-V(i)X)\ < a|0|iv0|[c|l- 

' C6'^\V(f>\\C\^\AC\^ (ifn = 2) 
C6^\V4>\\C\-4\AC\-4 (ifn = 3) 



< 



A7 



(4.65) 



|A((vV)C,AC)| < AIv|^4IVC|l4|AC| 

' CA|vli4lVC|^|AC| 



< 



< 



(if n = 2) 
C7A|v|i4|VC|^|AC|3 (ifn = 3) 

^|ACp + C!v|4,[VCl2, (ifn = 2) 
^|AC|2 + C|v|i,lVCp, (ifn = 3) 



(4.66) 



|A((0-V)0,AC)| < A|0|i4(p^)jV0|i4(p^)jAC|<A54|V0|^4(r,)|AC| 

At, 



< 



■|AC|^ + C752|V</.|i4(r,), 



(4.67) 



< 



^|AC||C|Le|'/'' + <^V' + V''lL3 



< ClACllClHid^lie + IVlie) 

< ciacIICIhi(I</'Ihi + IV'I?,i) 
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< ^\AC\' + C{\<P\%, + \ij\j,,)\C\l^. (4.68) 



Adding (|4.56p . (|4.57p and (|4.58p . then using the estimates above and taking 5 = cv, we 
finally obtain 

Ijt ('"^''^^^ + - + '^"^'^^ + C\Ai;\l) |wp + C (1 + |v|i, + 1^1? + \^P\i) Id? 

+ iy{C+ |Vvp + [Vvlio. + C\V^\^) + Cu^ (l + |vUoo(r,)) I Vu|i2(r,) 
+ C^.|Vu|i,(r,)+^^^^|V0li4(r,). (4.69) 

Denote 

y(t) = |w|2 + A|C|?, 

ait) = zy(C+|Vv|2 + |Vvlioo+C|V</.|2), 

6(t) = (1 + Iv|i^(r,)) |Vu|i2(r^) + Cz^lVul^^^^^^ + Cz/^IV./)!^^^^^). 

Since (v,-;/') is a smooth local solution on [0,T=k] and (p is bounded in L°°(0, T*; i7-'^(J7)), the 
inequality (|4.69p becomes 

<CY{t) + a{t) + b{t), (4.70) 
where a{t) < Cu and y(0) = [^^(O). Recah ()i:T7l) that |0|(t) < Cz^^s. Hence, 



sup y(t) < e^'^* (\e\^{d) + CvT^+ I b{t)dt] 

0<t<Tt \ Jo J 

< e^^* (^Cv + Ci^n + ^ * 6(t)dt^ , (4.71) 



where 



< CtJ ' \Vu\l2^^^)dt^ \cv ' \Vvi\l2f^^^^dt + Cv^, (4.72) 
where we used (12^61) and the fact H^{Vt) C L^{Q). Assume 

V / |Vu||2(r,)C^t ^ as ^ 0. (4.73) 

J 

Then combining (j4.7ip and (|4.72p . we have 

sup Y{t) ^ as 1/ ^ 0, (4.74) 
o<t<r* 

which equals to 

sup (|u — v| 2(t) + A|(/>-V|^i(t)) < sup y(t) + |0p(t) ^ as 1/ ^ 0. (4.75) 

0<t<T, 0<i<r, 
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In 2D case, if (u, (/>) is a global strong solution to ()l.ip - ()1.5p . then |Vup(t) is uniformly 
bounded in [0,oo), thus the assumption ()4.73p is automatically satisfied. (j4.72p is simplified to 

/ * bit)dt < for J, ^ (4_7g) 

Then, using (|4.76p . (|4.7ip becomes 

sup Y{t) < e^^* (Cu + Cun + Cv^l^\ < Cv^l'^ for v G (0, 1). 

0<i<T. ^ ' 

Therefore, 

sup (|u-v|2(t) + A|(/)-Vlii(t)) < sup + < C7zyi/2 for ^ g (0,1). (4.77) 

0<i<T. 0<t<T. 

The proof is complete. 



5. Related Models and Results 

Our system (|l.ip - ()1.3p is closely related to liq uid crystal model, Magnetohydrodynamics 



(MHD) equations, and viscoelastic system (see 28|] and the reference therein). In this section, 
we apply our previous approach to these models. 

5.1. Liquid Crystal 

If taking cj) as a vector, say d, system p.ip - ()1.3p can model the motion of liquid crystal flows: 

ut + (u-V)u = -Vp + i/Au- AV • (Vd(g) Vd), (5.78) 
V-u = 0, (5.79) 
dt + (u-V)d = 7(Ad-/(d)), (5.80) 

where d represents the director field of liquid crystal molecules. 

We consider the above system in 2D in an infinite channel M x (0, 1). We assume the velocity 
field u and director field d of liquid crystals are periodic with period 27r in the horizontal (x) 
direction. Set fl = (0,27r) x (0,1). u vanishes and the liquid crystal molecules align in a fixed 
direction d* at the boundary of the channel (i.e. at y = 0, 1), 

u = 0, d = d*, on (0,27r) X {0,1}. (5.81) 

(|5.78p - (|5.80p is subject to initial data 

u|t=o = uo(x,y), d\t=o = do{x,y), inn. (5.82) 

The above system is a sirnplified version of the Ericksen-Leslie model for the hydrodynamics 



im IS a sirnpli: 

f. BflflQ 



of liquid crystals (cf. [2, y, la, llO|] ) . Following exactly the same arguments as in [ll|] , one obtains 
the well-posedness results of system (|5.78p - (j5.82p . 
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The corresponding invisid liquid crystal model is 

u(°) + (u(°) • V)u(°) = -VP - AV • (Vd(°) O Vd(°)), (5.83) 

V-u(°) = 0, (5.84) 

+ (u(0) . V)d(o) = 7(AdW -/(d(0))), (5.85) 

u(°),d(°) are27r-periodic inx, u(°) • n = 0, d(°) = d*, on (0, 27r) x {0, 1}, (5.86) 

uW|t=o = uo(x,2/), d(o)|t=o = do(x,y), inJ7. (5.87) 

Following almost the same proof as that in Theorem 12.21 and Theorem 12.31 we have the follows. 

Theorem 5.1. Let s > 3 be an integer. Assume uq G Xs, do € H^^^. Then there exists a 
constant T,, > depending on \uq\s and |do|s+i, such that system (|5.83p - ()5.87p has a smooth 
solution (u^^), P, d^'')) in [0,T^] satisfying 

G L^{0,T,;Xs), P G L°°(0, T,; i7^+i(J7)), d^") G L^{0,n; H'+'). 

Moveover, (u'-''-*, d^'^)) is unique. 

Theorem 5.2. Let (u, d) be a global strong solution to (|5.78p - (|5.82p . and {u^^\d^'^^) a local 
smooth solution to (|5.83p - (|5.87p on Q x [0,T^,]. Then 

|u(t) - u(°)(t)|2 + A|d(t)-d(°)(t) I? ^0 as z^^O, VtG[0,r*]. (5.88) 

The convergence rate in ()5.88p is ci^"*^/^ for small v G (0, 1). 

5.2. The MHD equations 

The nondimensional form of viscous incompressible MHD equations is (cf. (isl ) 

ut + (u • V)u - S(B ■ V)B + Vp = ^Au, (5.89) 

ite 

Bi + (u- V)B- (B- V)u = TT-^B, (5.90) 

Km 

V-u = 0, V-B = 0, (5.91) 

where u, p and B denote the fluid velocity, its pressure and the magnetic field, respectively. Let 
L^,Uit:, B^, be the characteristic values for lengths, velocities, magnetic fields and densities of 
the fluid. The Reynolds number Rg = where u is the kinematic viscosity; the magnetic 

Reynolds number Rm = L^:U^,ap where /x is the magnetic permeability and a the resistivity of 
the fluid; S = p^i, ( = —^Yn] , where M is the Harmann number. 



Due to V • B = 0, we know in 2D case there exists a scalar function (j) such that 

Hence the above system in 2D is equivalent to the following equations 

uj + (u-V)u + Vp = ^Au-5V- (V</>(g)V0), (5.92) 

Re 
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V-u = 0, (5.93) 
0t + (u-V)</. = -^A0. (5.94) 

Compared with the system (|l.ip - ()1.3p . the nonhnear term /(0) vanishes in ()5.94p . 

Let 1/ — > and other parameters remain constants, we have Re ^ +oo and S, Rm are 
constants. System (j5.92p - (|5.94p formally becomes 

vt + (vV)v + VP = -SV- (VV'^ VV'), (5.95) 
V-v = 0, (5.96) 
V't + (v-V)V = ^A^. (5.97) 

We consider the vanishing viscosity limit of system (|5.92p - (|5.94p in the entire plane as 
i?e — > +0O. Omitting the terms with boundary layer function in the proof of Theorem [231 
and using the same estimates, we easily get similar result as follows. 

Theorem 5.3. Assume $7 = M^. Let (u, (^) he a global strong solution to ()5.92p - (|5.94p . and 
(v, ^/;) a local smooth solution to ()5.95p - ()5.97p on x [0, T*]. (u, 0) and (v,'i/;) are equipped 
with the same initial data. Then 

|u(t) -v(t)p + A|(/)(t) -V(t)|? ^ as i?e ^ +00, VtG[0,r*]. (5.98) 

The convergence rate in ()5.98p is cR~^ . 

Remark 5.1. In the literature, inviscid limit for MHD equations has been studied as both R^ 
and Rjn go to infinity, see 26, 27]. But in some industrial cases (see for example W]), Re ~ 
10^,i?m ~ 10^^,5" 1. Hence it is also meaningful to consider the inviscid limit for MHD 
equations as only Re — )• +oo. 

5.3. Viscoelasticity 

We consider the following system describing the motion of incompressible viscoelastic fluids 



(cf. 0,ll7|): 



ut + u-Vu + Vp = uAu + XV -{^^^^T^), (5.99) 
V-u = 0, (5.100) 
Tt + u-WF = VuJ". (5.101) 

Here, u represents the velocity field of materials, p the pressure, the kinetic viscosity constant, 
and A a parameter representing the competition between kinetic energy and elastic energy. 
is the deformation tensor, and W{J^) the elastic energy functional. 

Remark 5.2. The dynamic of viscoelastic fluids can be described by the flow map. Let X be the 
original labeling (Lagrangian) coordinate of the particle, and x the current (Eulerian) coordinate. 
For a given velocity u(x,t), the flow map (particle trajectory) x{X,t) is defined by the ODE : 

xt = u(x(X, t), t), x(X, 0) = X. 
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The deformation tensor P{X,t) is defined as P{X,t) = where we use the notation Tij = 
In the Eulerian coordinate, we can define T{x,t) = T{x{X,t),t) = T{X,t). The equa- 
tion ()5.10ip is deduced by the chain rule. The equation ()5.99p can be derived by the energetic 
variational approach, see fl^ . 0/ for intrinsic mechanics in detail. 

Taking divergence of both sides of (|5.1Uip and using V • u = 0, we get the transport equation 
for V • J-" as 

(V-7')t + u- V(V-^) = 0. (5.102) 
Since J-" is the identity matrix at t = 0, V • J^\t=o = 0, then (jS.lO'ip tells that 

V-7- = 0, yt>0. (5.103) 

In 2D, this implies there exists a vector function cf) = ((/>i,(/)2), such that 

If we take W{J-) as the Hookean elasticity energy, i.e., W{T) = \tr{FF'^), and use (ICTHD . 
then the system (|5.99p - (|5.10ip is equivalent to 

ut + (u-V)u = -V^ + z^Au- AV- (V(^«)V(^), (5.105) 
V-u = 0, (5.106) 
</)i + (u-V)0 = 0, (5.107) 

where the basic energy law is 

^|(|u|2 + A|V0|2) = -HVu|2. (5.108) 

Due to the system has partial dissipation, it is difficult to obtain global classical solution 
in general. Consider the system (j5.105p - (j5.107p in the entire plane R^. The existence of global 
smooth solutions near equilibrium has been proved in [3]. Moreover, the authors proved the 
local existence of a classical solution near equilibrium to the inviscid system in R'^: 

ut + (u-V)u = -Vg- AV- (V0 V0), 

V-u = 0, 
0i + (u-V)(/) = 0, 

by rewritting the above system as a quasilinear first order system. 

Due to the lack of dissipation term in the equation for (p, we cannot show the validity of 
vanishing viscosity limit for system (|5.105p - (|5.107p . by similar proof to that of Theorem 12.31 
But we believe that the inviscid limit holds for solutions near equilibrium. We are going to 
investigate this problem in forthcoming work. 
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